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We consider the free energy W[J] = Wt(H) of QCD coupled to an external source J^(x) = 
cos(£ • x), where is, by analogy with spin models, an external "magnetic" field with a 
color index that is modulated by a plane wave. We report an optimal bound on Wf,(H) and an 
exact asymptotic expression for Wt(H) at large H. They imply confinement of color in the sense 
that the free energy per unit volume Wk(H) /V and the average magnetization m(k,H) = i 
vanish in the limit of constant external field k — > 0. Recent lattice data indicate a gluon propagator 
D(k) which is non-zero, D(Q) ^ 0, at k = 0. This would imply a non-analyticity in Wk(H) at k = 0. 
We present a model that is consistent with the new results and exhibits (non)-analytic behavior. 
Direct numerical tests of the bounds are proposed. 
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1. Introduction 

Recent numerical studies on large lattices of the gluon propagator D(k) in Landau gauge in 
3 and 4 Euclidean dimensions, reviewed recently in [ffl], yield finite values for D(0) |0] - [f7|], 
in apparent disagreement with the theoretical expectation that D(0) = 0, originally obtained by 
Gribov [|p, and argued in The argument []9p which leads to D(0) = 0, relies on the hypothesis 
that the free energy W(J) in the presence of sources / is analytic in J at low momentum k. That 
hypothesis should perhaps be dropped in view of the apparent disagreement with the lattice data. 
This is of interest because a non-analyticity in the free energy is characteristic of a change of phase. 

The free energy W(J) enters the picture because it is the generating functional of the connected 
gluon correlators. In particular the gluon propagator is a second derivative of W (J) at J = 0, 



at , ^ 8 2 W(J) 



(l.i) 

7=0 



8J«{x)8J b v (y) 

The free energy W(J) in the presence of sources J is given by 

expW(J) = (exp(/,A)) 

= / dA p(A)exp(J,A), (1.2) 

JQ. 

where /x, V are Lorentz indices, and a,b are color indices, and 

(J,A) = j d D xjf l (x)A l l l (x). (1.3) 

The integral over A is effected in Landau gauge d^A^ = 0, and the domain of integration is re- 
stricted to the Gribov region Q., a region in A-space where the Faddeev-Popov operator is non- 
negative, M(A) = —d jl D jl (A) > 0. We use continuum notation and results, but we have in mind the 
limit of lattice QCD in the scaling region, that is gauge-fixed to the Landau (or Coulomb) gauge by 
a numerical algorithm that minimizes the Hilbert norm squared [|A|| 2 , and thereby fixes the gauge 
to the interior of the Gribov region. The vector potential, given by A{x) = gA pert (x), is unrenormal- 
ized, and has engineering dimension in mass units [A(x)] = 1 in all Euclidean dimension D, while 
[H] = D — 1 . (Our results also hold in the Coulomb gauge at fixed time, in which case D is the 
number of space dimensions.) The density p(A) is a positive, normalized probability distribution 
with support in the Gribov region £2. Because there are Gribov copies inside £2, p(A) is not unique 
and, in general, depends on the minimization algorithm. 
We consider a source that has the particular form 

/J(x) =#Jcos(fcti), (1.4) 

where we have aligned the 1-axis along k, so the free energy 

e\pW k (H) = ( exp[ J d D x ffJcos(fcci)Aj(x)] ), (1.5) 

depends only on the parameters k and H b . This is sufficient to generate the gluon propagator for 
momentum k, 

D°?(k)=2 d ' W f\, (1.6) 

1 J 
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where 

/ x W k (H) 

w k {H) = -*±-L (1.7) 

is the free energy per unit Euclidean volume. Because A»(jc) is transverse, only the transverse 
part of H is operative, and we impose k^H^ = 0, which yields H" = 0, and we write Hf, where 
i = 2,... D. By analogy with spin models, Hf may be interpreted as the strength of an external 
"magnetic" field, with a color index b, which is modulated by a plane wave cos(fcti). (This external 
magnetic field Hf, with color index b, should not be confused with the Yang-Mills color-magnetic 
field FPj.) 

A rigorous bound for W k (H) on a finite lattice was given in [^] which holds for any (numerical) 
gauge fixing with support inside the Gribov region £2. One can easily show that in the limit of large 
lattice volume V, and in the continuum limit, this implies the Lorentz-invariant continuum bound 
in D Euclidean dimensions, 

w k (H) < (2Dk 2 ) 1/2 \H\, (1.8) 



where \H\ 2 = Y.pL,b{^fi) 2 - A model satisfying the bound (L8) was recently exhibited in [ |i0| ] 



More recently, a stricter bound for w k (H) at finite H was obtained [11], that also holds for any 
(numerical) gauge fixing with support inside the Gribov region O, 

w k (H) < 2~ x l 2 k tr[{H a H a ) l l 2 }. (1.9) 

Here H a H a is the matrix with elements H"Hj. It has positive eigenvalues, and the positive square 
root is understood. A proof may be found in the Appendix that this bound is stricter than the old 



bound ( |l.8[ ). This bound is in fact optimal for a probability distribution p(A) of which it is known 
only that its support lies inside the Gribov region. 

The same expression also provides the asymptotic form of w k {H) at large H, and infinite 
Euclidean volume V flTij ] for any numerical gauge fixing with probability density p (A) with support 
that reaches all boundary points of CI, (but which may vanish on the boundary, p (A) = for A G d£l) 

w Kas (H) = 2- 1 / 2 ktr[(H a H a )^ 2 }. (1.10) 

Either bound yields in the zero-momentum limit 

w Q (H) = \imw k (H)=0. (1.11) 

As discussed in [g], this states that the system does not respond to a constant external color- 
magnetic field no matter how strong. It is a consequence of the proximity of the Gribov horizon 
in infrared directions. We shall return in the concluding section to the physical implications of this 
result for confinement of color. 

If wjc(H) were analytic in H in the limit k — > 0, eq. ( |1.11| ) would imply that all derivatives of 



the generating function wq{H) vanish, including in particular the gluon propagator (L6) at k = 0, 
D(0) = 0. However, as noted above, this disagrees with recent lattice data which indicate a finite 
value, D(0) ^ 0, in Euclidean dimensions 3 and 4. If this is true, then w k (H) must become non- 
analytic in H in the limit k — >■ 0. In order to get some insight about this, we examine the behavior 
of an improved model that has the exact asymptotic behavior ( |1.10| ). 
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2. Improved model 



The model is defined by the expression for the free energy 

r r 2 ua jja v 1 n -, r , k 2 H a H a \ 1 / 2 



w, mod (//)=^){tr[(/ + _ ? -) -/]-trln[2- 1 (/ + ™-) +2^/]}, (2.1) 



where > is an as yet undetermined function, and H a H" is the matrix with elements HfH", for 
i,j = 2,...D. This model possesses the following desirable features [pT|]: (i) It satisfies Wk, m od(ty = 
0, which is correct at H = for a normalized probability distribution JdA p(A) = 1. (ii) It has the 
asymptotic limit 

w kas (H) = lim W ^ mod(At//) = 2" 1 / 2 * tr[(M a ) 1/2 ], ( 2 - 2 ) 

that is correct at large H for any numerical gauge fixing that is strictly positive in the interior of the 
Gribov region £2. (iii) It satisfies the optimal bound 



Wk. 



mod (H)<2- l / 2 ktr[(H a H a ) 1 / 2 }, (2.3) 



which implies that the generating function vanishes at k = 0, wo, m od(#) = 0. (iv) The matrix of 
second derivatives is positive in the sense that 

d 2 w k d (H) h 
' dHfdH) 3 ~ 

holds for all v" and H", as required for this matrix to be a covariance. 

Because of the property wo^ mod (H) = 0, there must be some non-analyticity if, as indicated by 
numerical calculations, the gluon propagator D(k) at k = is positive D(0) > 0. It is instructive to 
see what kind of analyticity this would be in our model. Let A (H) > be the largest eigenvalue 



of the matrix HfHj. Inspection of (2T) shows that Wk, m od(H) is analytic in H inside a radius of 
convergence 

A(20 = ^. (2-5) 
Moreover from ( |2.1[ ) we have at small H, 

w k , mod (H) = J^H«H« + g(k)Olk 4 H 4 /g 4 (k)]. (2.6) 

For the gluon propagator D(k) ~ ^ sh^Th? ~ k 2 /§(k) to t> e finite at k = 0, as suggested by the 
lattice data, we must have g(k) = const k 2 near k = 0. In this case the coefficient of the H 4 term is 
of order l/k 2 , which diverges as k — > 0, as do all higher order coefficients. Moreover the radius of 
convergence of the series expansion of Wk,mod(H) is X(H) = 0(k 2 ), which vanishes like k 2 . 

Suppose that g(k) has a power law behavior g(k) ~ k v at k = 0. Then the radius of convergence 
behaves like X{H) ~ k 2v ~ 2 , which approaches with k for v > 1. The gluon propagator behaves 
like D(k) ~ k 2 ~ v , and w^ mo d(#) is non-analytic in H at k = when the propagator has a power 
law D(k) ~ with p < 1. Gribov's original calculation gave D(k) ~ k 2 /m 4 which corresponds 
to g{k) = 0(m 4 ), and w^mod^) is analytic in // at k = 0, with a radius of convergence = 
0(k- 2 ) — ^ oo for & — s- 0. 
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3. Conclusion 

By analogy with spin models, we define, for each momentum k, the analog of the bulk mag- 
netization in the presence of the external "magnetic" field Hf, 

*HM0 = ^, (3-D 

which describes the reaction of the spin system to the external color-magnetic field. Its physical 



meaning in gauge theory is apparent from QL5J ) which yields, 

M b (k,H) = {Jd D x cos(fcci)Af(x))fl 

= {\/2){a b {k) + a b {-k)) H . (3.2) 

Thus the "bulk magnetization" is in fact the k-th fourier component of the gauge field in the pres- 
ence of the external magnetic field. We also define the magnetization per unit (Euclidean) volume 

. Mf(k,H) 

mf(k,H) = lV v \ (3.3) 



given by 



<(*,*) = (3.4, 



The asymptotic free energy (1.10) determines the asymptotic magnetization per unit volume 
at large H, 

<„(*,") = %P 

= 2-^ 2 k[(H b H h )- l %Hj. (3.5) 

Its magnitude is given by (m" as m" af ,)(k,H) = k 2 /2, and we obtain the simple formula 

lim (mfmf)(k,H) = k 2 /2, (3.6) 

which holds for any numerical gauge fixing with support extending up to the boundary of the 
Gribov region Q.. 

We arrive at the remarkable conclusion that in the limit of constant external magnetic field, k — > 
0, the color magnetization per unit volume vanishes, no matter how strong the external magnetic 
field, 

lim lim m?(k,H) = 0. (3.7) 

Thus the system does not respond to a constant external color-magnetic field. In this precise sense 
the color degree of freedom m b (k,H) = ^(a b (k) +a b (—k))H is absent at k = 0. This conclusion 
holds whether or not the free energy Wk(H) is analytic in H in the limit k — > 0. Lattice data would 
indicate that it is not analytic. Besides reporting this result, we have presented a model, defined in 
(2.1), which saturates the asymptotic limit ( |1.10| ), and exhibits confinement of color. As we have 



seen, Wt, mo d(#) may be either analytic in H, or not, at k = 0, depending on the behavior of g(k) 
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at k = 0, but in either case, the conclusion stands, that the constant color degree of freedom of the 
gauge field is confined. 

Equations ( |L~^ ) and ( |1.10| ) may be checked numerically, at least in principle, by using the 
formula expWk(H) = (exp[f d D x //fcos(fcq)Af (x)]) to make a numerical determination of the 
generating function itself. For large values of H this may fluctuate too wildly. Alternatively one 
may measure the magnetization from the formula M^(k,H) = (J d D x cos(fcq)Af (x))#, where the 
source term Hj' cos(kx\)A^ (x) is included in the action that one simulates. This requires simulating 
the theory fixed in the Landau gauge instead of generating an ensemble from the gauge-invariant 
Wilson action then gauge fixing. It may be convenient to do this by numerical simulation of stochas- 



tic quantization [13] because that avoids calculating the Faddeev-Popov determinant explicitly. 
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A. The new bound is stricter than the old 



To prove that the new bound is stricter we write the old bound (|1.8|), for the case that k is 
aligned along the 1-axis, as Wk{H) <B\, where 

Bi = (2D) l l 2 k[\x(H a H a )] 1 / 2 . (A.l) 

We now diagonalize the matrix H"H" by a rotation in the (D — 1) -dimensional space, with eigen- 
values Xi, so 

Bi = (2D) 1 / 2 ^(£A ! ) 1/2 . (A.2) 

V i=2 

The new bound ( ]1.9| ) reads m>u{H) < B2 where 

B 2 = 2-^k^Xl'\ (A3) 

i=2 

and we must show B 2 <B\. (Note that the operations of square root and trace are interchanged.) 
We have 

D 

B\ = 2Dk 2 Y J X i , (A.4) 

i=2 

and 

B \ = 2- i e f a/ /2 a; /2 

<J=2 

\'=2 i^j=2 
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< 2- l k 2 



£ X i + 2- 1 £ (A« + A y ) 



(=2 1^=2 



D 



= 2-^ 2 (D-l) 



= (4D)" 1 (D-1)B 2 . 



(A.5) 



We thus obtain B\ < {AD) 1 (D — \)B\<B\, so the new bound is stricter, as asserted. 
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